Skyrmion Dynamics and NMR Line Shapes in QHE Ferromagnets 
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The low energy charged excitations in quantum Hall ferromagnets are topological defects in the 
spin orientation known as Skyrmions. Recent experimental studies on nuclear magnetic resonance 
spectral line shapes in quantum well heterostructures show a transition from a motionally narrowed 
to a broader "frozen" line shape as the temperature is lowered at fixed filling factor. We present 
a Skyrmion diffusion model that describes the experimental observations qualitatively and shows a 
time scale of ~ 50/^sec for the transport relaxation time of the Skyrmions. The transition is charac- 
terized by an intermediate time regime that we demonstrate is weakly sensitive to the dynamics of 
the charged spin texture excitations and the sub-band electronic wave functions within our model. 
We also show that the spectral line shape is very sensitive to the nuclear polarization profile along 
the z axis obtained through the optical pumping technique. 
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I. INTRODUCTION 

In the presence of a strong magnetic field, a two- 
dimensional interacting electron gas (2DEG) exhibits 
many different quantum states depending on field 
strength, electron density, and disorder. At Landau level 
filling factor i' — 1, when the number of electrons is equal 
to the number of available Landau orbitals in the low- 
est Landau level, the 2DEG is in an itinerant ferromag- 
netic state where all the electron spins are aligned with 
the magnetic field. El The novel features of this quantum 
Hall ferromagnet (QHF) state originate from the rela- 
tive strength of the electron-electron Coulomb energy, 
e^/eZs, and the Zeeman energy, g^sB, which are of or- 
der 160 K|-and 3 K respectively in typical experimental 
situations .□ The Pauli exclusion principle combined with 
the overwhelming cost of Coulomb excitations makes the 
single-particle spin-1/2 electron excitation gapped at fill- 
ing factor V — \. The low energy (but still gapped) 
charge excitations of the system are Skyrmion spin tex- 
tures containinffi,many (4-30) flipped spins and are topo- 
logically stable.El By paying the lower Zeeman price the 
spins can align locally creating a more advantageous 
charge distribution that lowers both the Hartree and ex- 
change Coulomb cost with respect to a single spin flip ex- 
citation. Because these excitations are the cheapest way 
to introduce charge into the system, they are present in 
the graU|nd state of the system at filling factors close to 
v = l.trH This produces a rapid reduction of the electron 
spin polarization away from u ~ 1 which is observed in 
experiments.l3 

The probe of choice to study the electron spin polariza- 
tion of the 2DEG is nuclear magnetic resonance (NMR) . 
Measurements of the Knight shift, which is linearly pro- 
portional to the electron spin polarization, have unam- 
biguously proven the existence of these exotic topologi- 



cal charge excitations.EI However, this rich probe keeps 
providing us with more surprising information about 
the Skyrmions, continuously challenging the theory of 
QHF.rJlecent experiments carried out by N. N. Kuzma 
et al.p using optical pumping techniques to enhance the 
NMR signal of the 2DEG, have focused on the NMR 
line shape as a function of temperature and filling factor 
near i/ = 1/3 (which is also a QHF state) . The observed 
free induction decay signals show a large dependence of 
the electron polarization on filling factor, indicating that 
charged spin excitations are more impoctant near v — l/i 
than previously expected from theoryU Near this filling 
factor, as the temperature is lowered, the Knight shift in- 
creases with decreasing temperature (although having a 
local minimum) until it reaches a saturation level at very 
low temperatures. The spectral line shape ranges from 
the motionally narrowed regime, where the polarization 
shows a sharp peak due to the_average polarization seen 
by each nucleus, to the frozerJj regime, where the peak 
is much broader due to the presence of frozen spin tex- 
tures in the ground state and the nuclei seeing different 
electron polarizations depending on their location. The 
interesting regime is the intermediate one where the dy- 
namical time scale of the Skyrmions is comparable to the 
inverse frequency of the Knight shift. This is the regime 
which contains the greatest wealth of information about 
Skyrmion dynamics. The situation is complicated by the 
fact that the Knight shift varies strongly with position 
in the z direction across the quantum well. Nuclei at the 
edges of the well see a lower electron density and hence 
a smaller Knight shift. As a result, these nuclei can still 
be in the motionally narrowed regime when nuclei in the 
center of the well are already close to the frozen regime. 
This subtlety is taken into account in our analysis. |-| 

Because the spin stiffness at = 1/3 is so smallQ the 
Skyrmions will be very small and the continuum field 
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theoretic approach will be poorly controlled. We there- 
fore focus on the case of filling factors g^f-r j/ = 1 where 
measurements are currently underway.B'cl We present a 
simple model for the Skyrmion dynamics near ^ = 1 in 
this paper. We find that the spectral line shape, besides 
being dependent on the transport relaxation time of the 
Skyrmions, is very sensitive to the nuclear spin polariza- 
tion density along the z direction and, to a lesser extent, 
the electron density profile in the z-direction. 

We organize this paper as follows. In Sec. II we in- 
troduce the theoretical background needed for our model 
calculations. In Sec. Ill we present our model and re- 
sults. In Sec. IV we discuss the implications of the results 
and possible new outlooks on this problem. 



II. THEORY 

In the free induction decay NMR experiments, af- 
ter the the wells are polarized by optical 



pumpingjO'Ej the nuclear spiriSpare tipped by a 7r/2 pulse 



and allowed to precess freely.a'cl These spins will precess 
at the Larmor frequency produced by the local magnetic 
field. This local magnetic field is composed of the exter- 
nal one plus a contribution from the electron polarization 
due to the Fermi contact hyperfine coupling which enters 
the nuclear spin Hanailtonian in the same way as the ex- 
ternal magnetic fieldEJ 



'^N = —gNi^N ^ • (Ho + Be(Rj)) 



(1) 



where [itq is the nuclear magneton, Sj is the nuclear spin 
in units of fi at position Rj, Ho is the applied mag- 
netic field, and Be(Rj) = (— 167r/iB/3) ^ ■ S^(5(ri — Rj) 
is the effective local magnetic field contribution due to 
the electronic polarization with and r.^ being the spin 
and position of the ith electron. 

The last term in eq. (Q) is the one responsible for the 
observed Knight shift. For the purposes of computing 
the Knight shift, it is adequate to replace Be by its ex- 
pectation value 



(s,(R,i/,r)) (X |w(R)|2p(R,z.,r) 



(2) 



where |m(R)P is the electron envelope function obtained 
from a self-consistent local spin-density approximation 
calculation, and ^(R, i^, T) is the average electron spin 
polarization at position R for a given filling factor and 
temperature. If we assume that the electron enve- 
lope function is only a function of z (the growth di- 
rection) we can further parameterize the local Knight 
shift as Ks{z) = pe{z)KsV{'R,i^,T), where Pe{z) = 
J dx dy\u{r)\'^ is the electron density along the z-direction 
(normalized to unity at its maximum) and Ks is a con- 
stant that can be fitted to the experimental spectra at 
the lowest temperatures. 



Next, we must connect these expressions to the ob- 
served intensity spectrum I {to), which is the time Fourier 
transform of the induced voltage produced in the tipping 
coil due to the precessing nuclear spins in the quantum 
wells. In the absence of in-plane spin nuclear decay, for a 
given nuclear spin, the time evolution of the spin's expec- 
tation value (relative to the evolution with zero Knight 
shift) is given by 

{S+{t))^{S,,+zS,y) 

= 5+(0)expH^ dTksPeizj)V{Rj,iy,T,T)] , 

where now we ha«e-|allowed the electronic polarization to 
vary with time.lliH13 The induced voltage is proportional 
to the time derivative of Re [(S'j''(i))] , which is approx- 
imately — Im [iuJo{Sj' (t))] whenever ujq ^ Kg, with ojq 
being the bare nuclear precession rate. However, for nu- 
clear spins in a solid the in-plane magnetization decays 
as a gaussian due to nuclear dipole-dipole interactions. 
This yields the final expression for the intensity 
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where p^ir) = (I],- <^(i"- 
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-Rj)S'2j) 



(5+(r,t)) 



(3) 



is the polarization den- 



sity of nuclear spins (not the number density of nuclei), 
and (7 — Au;/(2v21n2) with Aw being the full width 
half maximum (FWHM) of the unshifted NMR signal. 
In the samples used in the experiments Aw = 27r3.5kIIz, 
and therefore a — 9.34msec~^. Any model describ- 
ing the observed spectra has to contain realistic esti- 
mates of V{'Rj,i',T,t), Pei'r), and pN(r)- Of these three 
the first twfl|-Stffi_t|he--nnes that have been studied most 
extensively.Slpjia'y-r 
Kuzma et al. 



E3 In a previous model presented by 
to describe the transition near v — 1/3, 
a two polarization model for Viv^T^t) was used, and pe 
and pn were approximated by sinusoidal shapes. Here 
we attempt to improve upon this model by carefully ex- 
amining the approximations used for each quantity. 

At zero temperature V{'Rj,v), the electron polariza- 
tion, is exactly unity at = 1 as mentioned in section I. 
However, as the filling factor goes away from from v = 1, 
Skyrmions begin to appear in the ground state of the 
system. These charged spin texture excitations reduce 
the polarization locally and a realistic approximation for 
the shape of the Skyrmion is needed to obtain a reason- 
able V(Rj,h'). There HaMe-heeiujnany studies on the 
shape of the Skyrmions .BLnl3'L3~ll3 Most analytical ap- 
proaches describing the Skyrmion excitations have taken 
the route of effective field theories such as the modified 
0(3) non-linear sigma model (NL(t),Lj where a Zeeman 
term and a CoulombJnteraction term are inserted in the 
classical NLc model.ErQ'tj The magnetization profile of a 
Skyrmion obtained from this theory, in the limit of zero 
Zeeman energy, is given by 
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where A determines the size of the Skyrmion, and rj_ is 
the projection of r on the plane (denoted by r hence- 
forth). At larger distances away from the Skyrmion's 
center the Zeeman term dominates and the magnetiza- 
tion gees as 1 — mz(r) ~ e^'^'^^ /r, with n — AAIb near 
;/ = l.t3 This analytical approach has been very success- 
ful in explaining qualitatively the physics of the QHF, 
however it has not been able to predict quantitatively the 
shape of the Skyrmion excitations at experimentally ac- 
cessible parameters. In GaAs/Ala;Gai-2;As heterostruc- 
tures, at around 10 T, the number of ppin flips, per 
unit charge introduced near = 1 iso ~ 3 — 4. Sin- 
gle Skyrmion microscopic calculations using techniques 
such as Hartree-Fockta (HEX exact diagonalization, and 
variational wave functionsjlH have been more successful 
at obtaining quantitative agreement with experiments. 
These, however, are not as transparent as the classical 
model in describing the physics behind the excitations. 
Also, Brey et aZJl3 have performed HF calculations of the 
skyrme crystals formed at filling factors near v — 1 which 
most accurately.|describe the spin polarization observed 
in experiments. ElEa The failure of the NLcr model to pre- 
dict the correct Skyrmion shape for small K originates 
in the truncation of the gradient expansion. The HF cal- 
culations are in essence self-consistent mean field calcu- 
lations with the order of the gradient expansion taken to 
infinity and hence are more successful at predicting the 
small Skyrnugns that change shape on a much shorter 
length scaleO As described in the next section, we use 
a phenomenological form for ruz (f ) fitted to the Hartree 
Fock calculations. 

The electron density along the z-direction Pc{z), is ob- 
tained from the electronic spin-split energy levels of the 
GaAs heterostructures in the presence of a strong tilted 
field (as is used in the experiments). These charge distri- 
butions are calculated using a local spin-density approx- 
imation. The density profiles of the two spin states and 
two lowest levels, measured with respect to the valence 
band for a 30 — nm-wide GaAs single quantum well, are 
shown in Fig. |^. The|-parameters have been chosen to 
match the experimentsEI done near v — 1. 

The nuclear polarization profile estimate pjq, comes 
from considering the experimental set up. After a train 
of rf pulses which destroys the thermally induced nuclear 
polarization, the sample is radiated with circularly po- 
larized (7+ light tuned to the band gap in the well. This 
excites 100% polarized electrons to the lowest unoccupied 
level (spin down) and in the recombination process the 
nuclear system absorbs part of the angular momentum 
transferred to the sample. Hence, the optical pumping 
can produce an enhancement m the nuclear polarization 
by as much as a factor of 100.0 This optical pumping is 
necessary, at least at higher temperatures, to make the 
NMR signal visible. The problem is that, even though 
the initial excitation process is understood, the recom- 
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FIG. 1. Electronic density using the local spin-density ap- 
proximation for a sample with Btotai = 7.05T, a well width 
of 30nm, and a tilt angle of Ouit = 28.5°. The energies are 
measured from the top of the highest valence band. Note that 
the tilted magnetic field smears out the locations of the node 
positions for states in the first excited subband. 

bination process is much more complicated. In such a 
process, Skyrmions and anti-Skyrmions are being created 
constantly and this may have an effect on the nuclear po- 
larization profile which is not well understood at present. 

One may ask why at such low temperatures (1.6 K) 
and being so deep in the insulating phase of the quan- 
tum Hall state one can even observe motional narrowing 
of the NMR line shapes. We can answer this questionJia. 
reverse. Given the information from the experimentsJ3'cl 
what longitudinal resistivity can we infer? The onset 
of the frozen regime indicates that DK~^ ~ '^Ikyr ^^'^ 
dn/dp, ^ n/l^E^ with D being the diffusion constant, 
TT-skyr and n the density of Skyrmions and electrons re- 
spectively, and AE the disorder broadening of the Lan- 
dau level which can be estimated to be at least of the 
order of ~ lOK and is possibly much larger.Ej This in- 
formation can be inserted in the Einstein relation for the 
conductivity 

r, dn 9 1 Km , ^ 

^--^^d^-^"--A^ (4) 

3 X 10-i2f2-i 

— ' 

which gives a lower bound on the inverse conductivity 
^l^^xx ^ 300Gr2|i'— 1|. Hence, we see that the dynamics 
of the Skyrmions can appear to be fast on the NMR time 
scale even deep in the insulating regime. 

III. THE MODEL AND RESULTS 

When calculating the spectral intensity from eq. (3) 
it is useful to first think about the different time scales 
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in the problem and how relevant each one is in calculat- 
ing the NMR spectral line shape. At the experimental 
fields used (7.05 T) the bare nuclear precession rate is of 
order 100 MHz, the extra precession rate created by the 
electron spins (the Knight shift) is of order 20 KHz, the 
nuclear spin-lattice relaxation rate ranges from 4 mHz 
to 45 mHz, and the in-plane spin relaxation rate due ta 
the nuclear dipolar coupling is approximately 3.5 KHz.EI 
Hence, in the calculations that follow, we shall ignore 
the spin-lattice relaxation rate and place our zero of fre- 
quency at the bare NMR resonance. Furthermore, since 
we assume perfect cubic symmetry of the GaAs crystal, 
we omit any effects on the spectral calculation due to any 
nuclear quadruple splitting. 

Rather than doing a full microscopic calculation for 
'P(R.,i',T) in the presence of disorder, we take a more 
modest aim and focus on the transition where the spec- 
tral line shape goes from the frozen regime to the motion- 
ally narrowed regime. For — 1| ^ 1 the ground state of 
the 2DEG ai-T = is believed to be in a skyrme square 
lattice state.E3 In practice, the skyrme lattice-js_inelted 
in most of the accessible temperature ranges.E3c3 How- 
ever, at the low experimental temperatures considered 
here (T< 4 K), although the long length scale correla- 
tions vanish, we expect the short length scale correla- 
tions to contain crystal- like features. To model this we 
introduce a Skyrmion square lattice with the unit cell 
size given by the appropriate filling factor. We take the 
magnetization profile of the Skyrmion to be: 

mzir) = , 

+ A^e 

which, by choosing A and a appropriately, can i:^em- 
ble closely the HF calculations previously done.E3 For 
— 0.96 we used a = 5.6 and A^ — 2.6. Although in 
the HF calculations mz{0) ~ —0.6 due to the zero point 
fluctuations, we find that this has little effect on the spec- 
tral line shape, and we hence keep our simpler functional 
form in the model calculations. Also, in any experimen- 
tal situation, the part of the spectrum due to this small 
region of fully reversed spins tends to be weak. To model 
the time dependence of ^(R, T, t), instead of allowing 
the Skyrmions to undergo correlated thermally induced 
motion about their lattice points, we make the whole 
lattice move together to simplify the numerical calcula- 
tions. Hence, the lattice is only allowed to move collec- 
tively in a random walk by performing a jump of average 
distance I with a probability dt/rj in the time interval 
t to t + dt. The diffusion time across a unit cell of size 
L, Tdiff = {L'^/1'^)tj, is held constant for a given tem- 
perature and Tj and I are varied to test the sensitivity 
of the spectral line shape to the microscopic details of 
the dynamics. The motionally narrowed regime occurs 
when Tdiff ^ and each nucleus experiences the av- 
erage electron polarization in the sample. In the limit 
''diff ^ the Skyrmions are spatially frozen on the 

experimental time scale, with random motion being rare. 
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FIG. 2. NMR spectra (on the left) &X. v = \ for differ- 
ent nuclear polarization profiles (on the right). Profile (a) is 
taken from the lowest unoccupied spin-split energy levels (spin 
down) , profile (c) is a constant polarization density truncated 
at the edges, and (b) is a combination of (a) and (c). 

This frozen regime of the spectral line shape, observed 
at the lowest temperatures, contains information on the 
actual shape and static distribution of the Skyrmion exci- 
tations and the nuclear density profile. The intermediate 
regime, Tdis ~ , occurs when the FWHM of the spec- 
trum reaches its peak value. In this regime we find the 
greatest sensitivity of the line shape to the choice of rj 
for a fixed Tdiff. 

The adjustable parameter Kg, is fixed by fitting the 
peak frequency of the spectral line shape to the one ob- 
served experimentally in the frozen regime. The temper- 
ature is gauged by Tdiff' and calibrated by the onset of the 
motionally narrowed and frozen regimes observed in the 
experiments. We approximate the nuclear magnetization 
density by the electronic density from the lowest unoccu- 
pied level (El| in Fig. 1), since the induced polarization 
due to optical pumping is proportional to the local elefe 
tronic density of the electrons excited in such process.ll2l 
This gives a much more accurate profile jthan the ones 
approximated by simple sinusoidal shapes .□ One striking 
result of our model is the high sensitivity of the spectral 
line shape to the nuclear magnetization density. This sen- 
sitivity is illustrated in Fig. ^ where we show different 
spectral line shapes at = 1 for several nuclear polariza- 
tion profiles. This result strongly suggest that the details 
of the line shape will be difficult to understand from first 
principles without a better microscopic understanding of 
the optical pumping process and its effect on the nuclear 
polarization profile, together with other processes that 
may be affecting this profile. 

In Fig. H we show the spectral line shapes at different 
diffusion times (temperatures). The maximum FWHM 
is obtained between 40 and 50 fj,s. Also, note that Kg, 
which corresponds to the maximum in the spectra, de- 
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FIG. 3. NMR spectra at = 0.96 for different diffusion 
times. Here the parameters used are = 2.6 and ol = 5.6. 
Note that here the global electron polarization, V, is fixed, 
and therefore the peak location is no longer an accurate mea- 
sure of V. 

creases monotonically as the temperature (or rJTg) in- 
creases. It is important to note that at temperatures 
where the Skyrmion dynamics begin to 'freeze' in the 
NMR time scales, the peak of the spectrum (what is 
usually understood by the Knight shift) is no longer a 
good measure of the global average electron polarization, 
V = J cPrV{r,v,T,t — 0). Instead, in order to avoid 
dynamic effects, one should measure the first moment of 
the spectra when measuring since 



duj LjjI{lo) cx 



dzpN{z)pe{z) 



KsV , (6) 



where here we have duj I{oj) = 1 rather than /max = 
1 as used in Fig. 3 and 4. 

The sensitivity of the line shape to the choice of tj 
at Tdiff ~ 40 — 50/iS (corresponding to the maximum 
FWHM) is shown in Fig. ^. This sensitivity, although 
weak, is completely absent at other Tdiff' 's in the other 
regimes. 



IV. CONCLUSION 

Our model illustrates qualitatively the behavior of the 
NMR spectra as a function of temperature for T < 4 K. 
It accurately predicts a peak in the FWHM as a func- 
tion of temperature. It also shows a monotonic increase 
of the Knight shift with decreasing temperature, reach- 
ing a plateau at the lowest temperatures (largest Tdiff), 
however, it does not reproduce the local minimum in the 
Knight skift as a function of temperature observed near 

= l/3u Since the spin stiffness near v = 1/3 is very 
small, we do not expect this model to be valid near such 
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FIG. 4. NMR spectra for Tdifi — 40/isec and different tj's. 
Here u = 0.96, = 2.6, and a = 5.6. 



filling factor, where the Skyrmions are small and not well 
understood as is the case near v — I. 

We have also demonstrated that a full understanding of 
the NMR line shape must involve a better understanding 
of the nonequilibrium nuclear polarization profile. This 
profile is affected primarily by the optical pumping and 
possibly by other thermal relaxation processes. We also 
have shown that to measure the behaviour of the aver- 
age electron polarization V at these temperatures, one 
must measure the first moment of the intensity spectrum 
rather than its peak. We emphasize that our model does 
not attempt to calculate the highly quantum-mechanical 
motion of the Skyrmions. This motion may involve a 
semi-classical percolation in the case of heavy nonlocal- 
ized Skyrmions or variable range hopping in the case of 
highly localized Skyrmions in a random potential. Our 
phenomenological model does capture the appropriate 
time scales in the experiments which should be an es- 
sential result of a more sophisticated treatment of the 
problem. 
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